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ABSTRACT 

In the Large Area Crop Inventory Experiment 
a technique was devised using a vector transfor- 
mation of Landsat digital data to Indicate when 
vegetation is undergoing moisture stress. A 
relation was established between the remote- 
sensing-based criterion (the Green Index Number) 
and a ground-based criterion (Crop Moisture Index). 


1 . INTRODUCTION 

Meteorological Indexes which provide methods of computing the extent and 
severity of general drought have been developed. The most common of these are the 
Palmer and Foley drought indexes [1]. The Palmer drought index is the more 
Involved and the most satisfactory solution to the problem of combining precipita- 
tion and temperature data and producing a predictor variable [2]. During periods 
when a major drought is developing and spreading, it provides a means for routinely 
assessing the areal distribution of the various degrees of drought severity [3]. 

Palmer used a two-layer soil moisture model for estimating soil moisture defi- 
ciency on a month-by-month basis [^*]. These monthly variations are adequate for 
the study of drought; however, they do not reveal the shorter periods of soil mois- 
ture fluctuations which are Important during strategic phases of plant growth. The 
earlier drought study was modified by Palmer to calculate a weekly index of soil 
moisture related to crop water requirements [3]. This index uses the same two- 
layer models to compute a weekly hydrological balance. It is published in the 
United States as the Crop Moisture Index (CMI). Although the CMI provides an indi- 
cation of current plant-water relationships, it requires climatological analysis of 
a long record in order to derive the constants which define the moisture character- 
istics of the climate in the area of Interest. This paper describes a procedure 
using criteria based on remote sensing to detect and monitor crop moisture defi- 
ciencies without analyzing a long record of climatological data. 


2 . METHODS 

The procedure described in this paper was devised in an attempt to quantify 
the subjective judgment of the analyst-interpreter in deciding that a region is or 
is not drought affected. The data used are Landsat multlspectral scanner (MSS) 
values for Large Area Crop Inventory Experiment (LACIE [5]) sample segments through- 
out South Dakota, which were acquired during the 1975 and 1976 crop years (fig. 1), 
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and the CMI values reported weekly during these crop years by the U.S. Department of 
Agriculture (USDA) for each crop reporting district (CRD). Using the Ideas pre- 
sented by Kauth and Thomas [ 6 ], we devised a summary number, the Green Index Number 
(GIN), which was computed for each acquisition of each segment. Each segment was 
classified twice as drought affected or normal — first using only the GIN (remote- 
senslng-based criterion) and second using only the CMI (ground-based criterion). 

The results of these two classifications then were compared (table 1 ), and a con- 
tingency table of the GIN and CMI classifications (table 2 ) was prepared. 


The GIN Is defined as follows: First, the data in the segment acquisition 

are summarized by clustering using the Iterative Self-Organizing Clustering System 
(ISOCLS) algorithm [ 7 ] as Implemented on the Earth Resources Interactive Processing 
System (ERIPS [ 8 ]) on the special purpose processor. (This parallel processor clus- 
ters a segment In approximately 30 seconds.) The clustering procedure summarizes the 
segment In 20 or fewer cluster means In the four Landsat channels. The count of 
pixels belonging to each cluster Is also calculated. Each mean vector x^ Is then 
transformed by 

y^ = Ax^ + b* ( 1 ) 

where 


y 


1 


A 


b' 



a vector representing the Landsat -1 version of the Kauth-Thomas trans- 
formation of x^ [ 6 ]; the subscript number indicates the Landsat channel, 
and the supercrlpt Is the cluster number. 


0.44326 

0.6325 

0.5857 

0.26441 

-.2897 

-.5620 

.5995 

.44907 

-. 82442 

.5329 

-.0502 

.1850 

.2229 

.0125 

-.54431 

.80944 


( 0 . 145 ,- 1 . 50 , 10 . 61 , 2 . 22 ) 


Each vector is inspected automatically, and any vector having values unreasonable 
for agricultural data is discarded using the following procedure. 


A cluster y^ is accepted as good only If 

30 < yj < no' 
-10 < y\ 

-10 < y3 
-10 < yj < 10 


( 2 ) 


The greenness level m of the soil line then Is estimated by the minimum second 
channel value y| for acceptable clusters. That is, 

ra = M I N (yj\ (3) 

i is good ' ' 

Then the green number g^ is computed for each cluster by 

g^ = y^ - m (44) 
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The value of g should be a good measure of the green vegetation present on the 
picture elements (pixels) In cluster 1. Experience with spectral plots In the 
brightness-greenness plane and the corresponding Imagery led to the following 
assumptions . 

• g^ = 0 Indicates bare soil 

• g^ = 5 Indicates a trace of vegetation 

• g^ = 15 Indicates good cover of vegetation 

Operating on these assumptions, we chose the level 15 and defined GIN to be 
the percentage of pixels In the entire Image within clusters having green numbers 
greater than 15. This value was considered somewhat unstable for conseoutlve-day 
data where, for example, a cluster would slip from 15.1 one day to 1^).9 the next 
day; therefore, cubic weighting was added to smooth this calculation In the follow- 
ing manner. 

Let a cluster have green number g^ and contain n pixels. Define the weighting 


factor w by 

w = 

0 If 

g^ 

< 11 





w = 

1 If 


> 17 


► 

(5) 

Otherwise, 

w = 

1 *^ 

1 






The cluster is counted as having w x n pixels with green numbers greater than 15. 
This curve makes a smooth transition from full counting to not counting as the 
green number decreases. 


The GIN then Is an estimate of the percentage of pixels in a Landsat scene 
having green numbers high enough (> 15 ) to indicate full cover of green vegetation. 
It Is computed using only Landsat data. A sample spectral plot of green numbers 
versus brightness [6] Is given In figure 2. 

Having defined this segment summary number, the results had to be analyzed. 

For the test set used. It was impossible to make (1) year-to-year comparisons 
because of sparse Landsat acquisitions and Important weather differences between 
1975 and 1976, (2) area-to-area comparisons because of differing wheat propor- 
tions and haze effects, and (3) a direct comparison of GIN's and CMI's because 
the GIN contains crop stage effects not taken into account In calculating the 
CMI. This comparison also Is confused by nonwheat crops with different crop 
calendars from those of wheat. 

Therefore, the preferred technique for this analysis was a comparison of clas- 
sifications. The plot of GIN versus time for a normal, predominantly wheat segment 
should follow a curve such as a In figure 3. If an observed point for a segment 
fell into the shaded region, the segment was classified as drought affected. The 
bounds for the shaded region were defined empirically as shown In figure 3, with t 
defined as the approximate spring emergence date In days. For different areas or 
years, the shaded area can be moved from side to side to match the greenup curve. 
The Initial point In South Dakota was usually near day 110 (t = 110). This 
classification was compared to a classification based on the CMI for CRD's, wherein 
a CRD was classified as drought affected If the CMI fell below -0.5 for 2 consecu- 
tive weeks. Both classifications were restricted to similar time frames. Classi- 
fication was performed only for data between April 1 and July 10. 


3. DATA AND RESULTS 

The data used In this study consisted of all LACIE segments In South Dakota 
which had at least 5 percent wheat as measured by the LACIE Classification and 
Mensuration Subsystem (CAMS [9]) in the 1976 growing season. This definition 
yields 17 segments with 34 possible classifications. Of the 34, 12 had either 
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Insufficient data during the growing season or data were inaccessible for other 
reasons. The final data set contained 22 segment yearsl for 13 LACIE segments 
(fig. 1, table 1). The contingency table (table 2), which applies the two classi- 
fication methods to the 22 good segment years, shows that the classifications based 
on the CMI and GIN are related. It was concluded that the GIN is detecting moisture 
through crop responses. 

An inspection of the five disagreements on the classification results (table 1) 
disclosed that on two segments the GIN algorithm was confused by a lake. Three 
of the segments were on the edge of their CRD's, and the CMI classifications may not 
reflect the actual conditions in these segments. One segment was located on the 
eastern edge of the CRD where heavy rains occurred at the weather station located 
in the Black Hills in the western part of the CRD, causing a possible Incorrect 
condition to be reflected by the CMI. 

Examples of the segment classification procedure are shown in figures 4 and 5- 
The GIN indicates that 1975 was normal for the entire crop season for segment J 
(fig. ^). In 1976, the GIN indicated that by May 2k there was moisture stress in 
segment J. This Indicates that the GIN detected vegetation moisture stress at the 
saime time as the CMI. Segment L (fig. 5) experienced drought as indicated by both 
the GIN and the CMI during the 1975 crop year. In 1976, the GIN indicated moisture 
stress on May 26, which was confirmed by the CMI. 


4 . DISCUSSION 

The results reported cover a first attempt at measuring drought conditions 
using Landsat data. Preliminary checks indicate that this procedure works in the 
U.S. Great Plains. More recently, we have developed a program using raw image data 
and avoiding the clustering procedure, enabling the generation of GIN's with any 
threshold. This new procedure will be utilized in continuing investigations. 


5 . CONCLUSIONS 

A technique was developed using Landsat digital data from 9-by-ll kilometer 
sample segments, which Indicates when agricultural vegetation is undergoing mois- 
ture stress. A relation between this technique, which utilizes remote sensing, and 
a ground-based criterion (the CMI) has been shown. Indications are that this pro- 
cedure may be capable of detecting crop moisture deficiencies in areas of the world 
where ground information is not available or reliable. 
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TABLE 1. RESULTS OP GIN AND 
CMI CLASSIFICATIONS 


TABLE 2. CONTINGENCY TABLE OP GIN AND 
CMI CLASSIFICATION METHODS 


D = Drought conditioni 
W = Normal conditions 
— = No data 


Segment 

1975 

1976 

GIN 

CMI 

GIN 

CMI 

A 

w 

W 

- 

D 

B 

- 

W 

W 

D 

C 

W 

W 

D 

D 

D 

W 

W 

D 

0 

E 

w 

w 

D 

D 

F 

w 

w 

W 

D 

G 

w 

w 

W 

D 

H 

w 

w 

D 

W 

1 

w 

w 

- 

W 

J 

w 

w 

D 

D 

K 

w 

w 

D 

D 

L 

D 

D 

D 

D 

M 

- 

D 

W 

D 


CMI 



Normal 

Dry 


Normal 

10 

4 

14 

Dry 

1 

7 

8 


11 

11 

22 


= 7.07 with 1 degree of freedom. 
P = 0.0082 ° level of si^iificartca. 


NORTH DAKOTA 



FIGURE 1. SEGMENT LOCATIONS. Map of South Dakota showing locations of 
LACIE 9-by-ll kilometer sample segments. 
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FIGURE 2. SPECTRAL PLOT. Sample spectral plot of cluster statistics. 



FIGURE 3* PLOT. GIN versus time for a normal, predominantly wheat segment. 


PERCENTAGE OF PIXELS WITH GREEN NUMBER ^ IS 
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